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Motivated in part by recent experiments on liquid crystals with bent-core molecules, which are 
observed to display a spontaneous chiral symmetry breaking, we introduce a field theory of a 3rd- 
rank tensor order parameter T'-''' to describe the isotropic-to-tetrahedratic phase transition that 
we predict to take place in these materials. We study the critical properties of the corresponding 
phase transition and find that this transition, continuous at the mean-field level, is generically driven 
Ist-order by thermal fluctuations. 



Experimental studies of liquid crystals, systems that 
exhibit phases that are intermediate (both in their prop- 
erties and the symmetries that they break) between ordi- 
nary liquids and solids continue to present new surprises 
and theoretical challenges. The experimental discovery 
in achiral bent-core molecules of polar-ordered smectic 
layers by the Tokyo Tech group and of spontaneous 
chiral symmetry breaking by the Boulder group |^ is the 
latest important example of such a challenge. Bent-core 
molecules are V-shaped with C2v symmetry Q charac- 
terized by a non-polar direction n pointing from one end- 
point of the V to the other and an orthogonal polar di- 
rection p pointing to the vertex of the V. Seven distinct 
phases of bent-core molecules tentatively labeled Bi to 
By have been identified Q|, though not all have been 
fully characterized. Two phases, the B2 and By phases, 
are smectic phases consisting of stacks of fluid layers with 
some internal tilt order. The Boulder group ||^ has shown 
that the B2 phase is an anti-ferrolectric smectic-C phase 
in which the configuration of molecules in each layer is 
chiral with chirality alternating in sign from one layer 
to the next. A material of achiral nematogens having 
a ground state which is ferroelectric and homogenously 
chiral has also recently been discovered. |Q Within smec- 
tic phases, such spontaneous breaking of chiral symme- 
try is possible because layers provide an anisotopic envi- 
ronment, defined by the layer normal, relative to which 
molecules can tilt and rotate. 

Though most bent-core molecules do not exhibit ne- 
matic phases, some do and it is of some interest to de- 
velop the minimal theory capable of exhibiting isotropic, 
nematic, and the various smectic phases of bent-core 
molecules. Homogeneous (i.e., translationally- invariant) 
phases composed of achiral molecules and characterized 
only by two orthogonal axes cannot be chiral. Clearly, 
then, a theory, whose orientational order is based only 
on the Maier-Saupe symmetric-traceless 2nd-rank tensor 



order parameter Qij cannot describe the exotic orienta- 
tional order of the B2 phase. In many liquid crystals, 
electric dipoles play an insignificant role in determining 
the development of order, and it is appropriate to con- 
sider models in which molecular properties are described 
entirely by their mass-moment tensors. Since the first 
mass-moment tensor relative to the center of mass is zero, 
these molecules do not have mass-moment vector order 
parameter. They do, however, have a third-rank mass- 
moment tensor that does contain information about the 
polar shape of a bent-core molecule. The simplest model 
capable of describing both the nematic and the orienta- 
tional order of the aforementioned smectic phases of bent- 
core molecules is one with both 2nd-rank and 3rd-rank 
order parameters. Arbitrary symmetric third rank ten- 
sors can be decomposed into a symmetric-traceless part 
T'^*^, which transforms under the I = 3 representation of 
the group 5'0(3), and a part that transforms like a vector 
p' {I = 1) under SO{3). 

A model with Q'-' , T*-''' and order parameters pro- 
duces an incredibly rich phase diagram at the mean-field 
level. 1^ It predicts, in particular, translationally invari- 
ant nematic phases with D^h, D2d and even chiral D2 
symmetry in addition to the usual uniaxial and biaxial 
nematic phases with respective Daoh and D2h symmetry. 
The latter phases are characterized by = 0, T^^^ = 0, 
and Q*-' 7^ 0, whereas the former, more unfamiliar phases 
have both T'^'^ and Q^^ nonzero, but p* = 0. There is 
also a novel tetrahedratic (T) phase with p' = 0, Q'^ = 0, 
and only T^i^ nonzero. The direct transition from the 
isotropic phase to the tetrahedratic phase can be de- 
scribed by a simpler model depending only on the 3rd- 
rank tensor order parameter T^i^ , Since there are no odd 
order invariants of T^i^ ^ this transition can be continuous 
in mean-field theory. ||^ In this note, we study fluctuation 
corrections to the mean-field isotropic-to-tetrahedratic 
(IT) transition within an e-expansion about the upper 
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critical dimension dc = 4. Our principal result is that 
the presence of two distinct fourth order invariants (see 
Eq. (||) below) leads to a runaway within the e-expansion 
whose effect in analogy with charged superconductors 
and a ferromagnet in a cubic field is to convert this 
mean-field 2nd-order transition to a Ist-order one. 

Our starting point is the Landau free-energy density 
functional T-l[Tijk], controlling the direct transition into 
the tetrahedratic phase: 

In the above expressions we have used the Einstein sum- 
mation convention, and have left out the dipolar-like gra- 
dient terms dij^Ti-^j^-di^Ti^jk that couple internal indices 
of Tijk to that of spatial coordinate x. This latter sim- 
plification, quite commonly employed in the analysis of 
the isotropic-nematic (IN) transition, might in principle 
modify the nature of the IT transitions at the longest 
scales, as it does for example in ferromagnets. Q We 
leave the analysis of these "dipolar" effects to a future 
investigation. We have also omitted couplings of T*-''^ 
to other order parameters. Since here we are interested 
in the direct transition to the T phase, all other order 
parameters (e.g., the nematic Qij) are generically non- 
critical and can, therefore, be safely integrated out, only 
finitely renormalizing the Landau parameters in 7i. As 
usual the quadratic parameter r T — Tc vanishes at 
the mean- field transition temperatures Tc, which is de- 
termined predominantly by the interaction potential be- 
tween molecules. In writing down quartic nonlinearities 
in 7i we have used a non-obvious relation. 
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(2) 



valid when the dimension dx of the space in which T'*J'= 
rotates is three, to reduce the number of independent 
quartic couplings in Ti.[T] from three to two. The cou- 
pling V was incorrectly omitted in the analysis of Fel. 
§ If w = 0, the model described by Eq. @) has 0(7) 
symmetry. Thus, it is analogous to an 0{N) Heisenberg 
ferromagnet with an A^-dimensional vector order param- 
eter Tq, in a cubic crystal field, giving rise to a v'^^T^ 
quartic potential in addition to the 0{N) invariant po- 

tentialu(E„7^a)'- 

As discussed extensively in Ref. j?], the IT transition 
is generically continuous at the mean-field level since ro- 
tational invariance forbids the appearance of a cubic in- 
variant in T^J'^, Its nature, however, depends on the sign 
of V. For V > 0, the stable low-temperature phase has 
Djh symmetry, whereas for u < 0, it has tetrahedral (T^) 
symmetry. Similar behavior is seen in 0{N) ferromag- 
nets in a cubic crystal field. When the cubic coupling v is 



nonzero, the low-temperature phase has Ta ^ (1, 0, 0) 
for u < and Tq, ^ (1, 1, 1)/^/N for u > 0. Our main 
goal then is to go beyond this mean-field theory [Q and 
study effects of thermal fiuctuations on this heretofore 
unexplored IT transition. 

Before analyzing fluctuations in this model, it is useful 
to introduce a tensor representation for the 7 indepen- 
dent components of T'^i^ _ Any third-rank symmetric- 
traceless field r*^'^(x) can be expressed as T^^'^{x) = 
ELl^^a(x)/^^^ where 
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are an orthonormal set of basis tensors with n% m*, and 
r forming a space-fixed right-handed orthonormal basis. 
The 0(7) invariance of the Eq. d) when u = follows 
easily from the relation T^if^T^J^ = Y.''a=i ^a- 

The effective Hamiltonian H ^ J d'^xTC in d spatial 
dimensions can be conveniently rewritten in terms of 
Fourier transformed T-fields 
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V;'^;^" is a fully symmetrized quar- 
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tic potential (consisting of products of 3 Kronecker 5- 
functions), containing the u and v vertex couplings and 
is too complicated to be reproduced here. There are at 
least two ways to generalize our model, originally defined 
in 3 dimensions to arbitrary spatial dimension d. We can 
either fix the dimension dr of the space in which T*-''^ 
"lives" to be three and let the dimension d of space be 
arbitrary, or we could lock dx and d together. Here, we 
adopt the former option, keeping dx — 3. If we had in- 
cluded dipolar- like gradient coupling terms in Eq. (Q), 
we would have had no choice but to use the analytical 
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continuation in which dx = d. The advantage of sym- 
metrized representation in Eq. is that averages of any 
pair of T*-''^ fields appearing in the quartic operator are 
equivalent, and the field theory part of the calculation is 
effectively mapped onto a standard analysis of a scalar 
(j)^ field theory. 

It is not difficult to show that a corresponding tensor 
part of the T — T propagator G'^" defined by 



(T,,fe(k)Ti™„(k')) = (27r)'^5('^)(k + k')^ 
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is given by 
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The tensor G-^" is uniquely constrained by the rota- 
tional invariance to be a sum of products of Kronecker 
(5-functions, with two nontrivial coefficients uniquely de- 
termined by the various traces of G'^". Another more 
systematic way of obtaining G-^" is to use the 0(7) rep- 
resentation of Tyfe(k), Eq. inside the left hand side of 
Eq. to produce 



(Ty7c(k)T/„„(k')) = 



/^■'=/iT"(r,(k)r,,(k')) . (7) 



a.Q' — 1 

Using the 0(7) representation of H along with equipar- 
tition, we obtain (Ta(k)Ta' (k')) = {271^ 5^'^\\<i + 
'^')5aa' / {k"^ + r), which when used inside Eq. (^ gives 
G'^" = ELi I'J^liT'" that can be shown to be equiva- 
lent to our expression for G-^" in Eq. (||). 

To assess the role of thermal fluctuations, we need to 
compute the total free energy by integrating over the 
field T'^^^{'x) rather than by simply minimizing the effec- 
tive Hamiltonian above, as we did in Rcf. 0. As usual, 
because of the nonlinear quartic interaction this cannot 
be done exactly, and near a critical point (r « 0), a 
perturbative calculation in V^^'^^ diverges (for d < 4), 
thereby indicating that the nature of the IT transition is 
qualitatively modified by thermal fluctuations. 

To make sense of these divergences we employ the 
standard perturbative momentum-shell RG calculation 
[^2|, in which we integrate out an infinitesimal frac- 
tion (a momentum shell) of degrees of freedom at a 
time near the shortest length scales, at the lattice cutoff 
A~^. More concretely, we separate the T^-fe order pa- 
rameter field into high and low-wavevector components 
Tyfe(x) = T^^fe(x) + T^^fe(x), with T>.^^ having support in 
the wavevector range Ae~^ < A: < A and A is an ultravi- 
olet (uv) cutoff of order of the inverse of molecular size. 



We integrate out the high-wavevector part T^^j,(x), per- 
turbatively in Vj'^:^*', and rescale the lengths and long- 
wavelength part of the fields according to x = e^x', 
T;<j,(x) = e^^^~'^'>/^T,jk{^'), so as to restore the uv cutoff 
back to A and to (1-loop order) keep the coefficient of 
the fc^|T!yfc(fc)p in H fixed at 1. Under this transforma- 
tion the resulting effective free energy functional H can 
be restored into its original form Eq. (^) with effective 
^-dependent couplings. 

3^4 
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FIG. 1. Graphical representation of the leading order cor- 
rection to the reduced temperature r and the quartic coupling 
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Standard perturbative 1-loop RG analysis, represented 
graphically in Figjl] leads to the following flow equations 



^G-i 



2G: 



12A'*-2Gd 
(l + r/A2) 



di '1''=^ 



"(l + r/A2)^ 



^^i'* Gi., i' Gi^ i' 



(8) 



(9) 



where e = 4— d. These equations describe the evolution of 
the reduced temperature r{£) (determined by the q — 
part of the 2-point inverse propagator Gr\^((j' = 0,r)) 

and the quartic vertex Vi^i^{£) after a fraction of 
high-wavevector modes has been integrated out. In above 
Cd = ((27r)''r(d/2)) is a surface area of a d- 

dimensional sphere divided by (27r)'', an obvious exten- 
sion of Einstein convention was employed, and e was 
taken to be small so as to allow higher-order corrections 
in V";'^;'^'' to be neglected. We note in passing that for the 

fully symmetrized y^'l^ vertex, the factors 12 and 36 of 
the respective second term in Eqs. ^ and (||) are identi- 
cal to the corresponding coefRcients in the equations for 
the reduced temperature and quartic coupling constant 
in the scalar 0"* (Ising) model. 

After technically tedious but conceptually straightfor- 
ward tensor contraction in Eq. (^), we find that the 
RG flow equations at criticality [r — 0) for the two di- 
mensionless couplings of our theory u = ut^^^'^Cd and 
V = vA^~'^Cd, with u and v defined in Eq. (m, are 
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The flow off the critical surface, controlled by f{t) 
r{i)/A'^ is determined by the equation 



—77^ = 2r + r 36u + —v) 

M 1 + r ^ 5 ' 



(11) 



As a simple che ck o n our results, we first note that the 
flow for w, Eq. (lOt) has an obvious fixed point solution 
of u = 0. This property of the flow for u is a rigorous 
consequence of an exact 0(7) invariance of our theory for 
■D = 0, as can be easily seen from Eq. (|l|). 

Standard fixed point analysis of the above equations 
leads to two physical critical points given by: (i) Gaus- 
sian with u = 0, w = 0, and (ii) 0(7) with u — e/60, -5 = 0. 
As a further check, we also note that for w = 0, above flow 



equations, Eqs. (gTj) and (10a) reduce exactly to the flow 
equations for the 0{N) Heisenberg model with N = 7 
(the model that our Hamiltonian H explicitly reduces to 
for V — 0), with corresponding universal coefficients 36 
and 60, in above equations, identical to the well-known 
4(iV+2) and 4(iV-H8) coefficient for the 0{N) Heisenberg 
model with N — 7. 

As expected and discussed above, for d < 4 the Gaus- 
sian critical point is unstable. For the special v = 
0(7) invariant subspace, this instability terminates at 
the well-known N — 7 Heisenberg model critical point 
characterized by (1-loop) exponents: = 1/(2 — 3e/5), 
^ — 2v and 77 = 0. 




FIG. 2. RG flow characterizing the IT phase transition. 
The Gaussian and 0{7) (it = e/60, u = 0) critical points are 
indicated by black circles. In a generic model, with d 7^ 0, 
the flow runs away into region of u < 0, suggesting a fluctua- 
tion-driven Ist-order IT transition. 

This 0(7) critical point is also unstable with respect 
to turning on v. Standard linear analysis of Eqs. ( p^ ) 
around the 0(7) {v = 0) critical point leads to two eigen- 
vectors and eigenvalues that characterize the instability 



to turning on the 0(7)-symmetry breaking perturbation 

V 

{5u, Sv)u = u{l , 0) , Af, = -e , (12a) 

(Su, Sv%^v{--^, 1), A, = (12b) 
90 5 

As is graphically displayed in Fig.|| for the critical sur- 
face rc{i), we find that in the full model the effective 
couplings u{i) and v{£) run off to infinity. 

The absence of a stable fixed point and the runaway to- 
wards negative values of u is analogous to the situation 
encountered in model superconductors ^ with {N/2)- 
component 17(1) symmetry in which there is runaway 
from an 0(A'^)-symmetric Heisenberg fixed point towards 
negative u when charge is turned on. In the latter sys- 
tem, the runaway can be identified with a fluctuation in- 
duced first-order phase transition PjT^]. It is likely that 
the runaway in the current problem is also a signal of a 
fiuctuation-induced first-order IT transition. In the case 
of the superconductor, an effective free energy indicating 
a first-order transition could be obtained by integrating 
out the vector potential. In the model considered here, 
there is no obvious analog of the vector potential. As a 
result, the detailed analysis of the first-order IT transi- 
tion is likely to be more similar to that of the fluctuation 
induced first-order transition, which exists for appropri- 
ate values of the potentials, in 0{N) magnets in a cubic 
field (l|], analyzed for TV = 2 by Rudnick 
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